double integration method for beam deflections